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Abstract: We investigated the floristic composition, woody species 
diversity and spatial distribution of trees in a tropical wet evergreen 
forest in Kaptai national park in Chittagong Hill Tracks, Bangladesh. We 
recorded 25 families, 37 genera, 40 species and 1771 woody individuals 
in a 0.09 ha plot. Euphorbiaceae and Moraceae were the most species- 
rich families, and Castanopsis, Ficus and Terminalia were the most 
species-rich genera. Bur sera s errata Wall, ex Colebr. was the dominant 
species in terms of highest importance value (13%). Trema orientalis (L.) 

B1 was typically a light demanding species as it appeared in the top can¬ 
opy with only one individual having the seventh highest IV, but had no 
regeneration. The expected maximum number of species (S max ) was 140, 
indicating that many species may invade the forest as the S max is greater 
than the recorded total number of species. The nature of the disappear¬ 
ance and appearance of species in the present forest reflects instability of 
floristic composition. The values of Shannon’s index H' and Pielou’s 
index j' (evenness) were 3.36 bit and 0.63, respectively. These values 
show moderately high species diversity as compared to other subtropical 
forests in the tropics. In addition, a sample area of 200 m 2 in this forest 
would be sufficient for measuring the diversity indices H' and J' , 
whereas the trend of J may indicate the rate of equality of individuals 
among the different species decreased with increasing area. The distribu¬ 
tion pattern for the total stand was completely random. However, the 
dominant species showed aggregate distribution for small areas, but 
random distribution for large areas. The spatial association between 
species showed that the strongest positive interspecific association oc¬ 
curred between Streblus asper Lour, and Castanea indica Roxb. (co = 
0.51). As a whole, most species were weakly associated with each other, 
of which 58% species associations were completely negative. The result 
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of cluster analysis showed that species pairs were spatially independent 
at all or most small clusters; stands of species from all clusters are mosa¬ 
ics of complete habitat and pioneer habitat. All patches in this forest 
community have similar habitat and regeneration niches, which could be 
a phenomenon for a young growth forest. Thus, the existence of habitat 
and regeneration niches may be an important factor in the maintenance of 
diversity in this forest. 

Key words: Chittagong Hill Tracts, floristic composition, spatial distri¬ 
bution of trees, tropical wet evergreen forest, woody species diversity 

Introduction 

Biodiversity is a key issue in nature conservation (Wilson 1988; 
1992) and species diversity is one of the important components 
of biodiversity (Ito 1997). Diversity of trees is fundamental to 
total forest biodiversity, because trees provide resources and 
habitats for almost all other forest species (Hall and Swaine 1976; 
Huston 1994; Whitmore 1998; Huang et al. 2003). Generally, 
tree species diversity in a forest varies greatly from place to 
place mainly due to variations in biogeography, habitat and dis¬ 
turbance (Whitmore 1998). 

Measures of species diversity play a central role in ecology 
and conservation biology. The most commonly employed meas¬ 
ures of species diversity are the Shannon function, species rich¬ 
ness (number of species), and evenness (the distribution of abun¬ 
dance among the species, sometimes known as equitability). In 
addition, spatial distribution of trees has been a major source of 
interest for plant ecologists because of its potential role in ex¬ 
plaining the coexistence of tree species in species-rich forests 
(Bunyavejchewin et al. 2003). Analyzing the spatial distribution 
patterns of trees may help to detennine the mechanisms impor¬ 
tant in structuring forest communities because such patterns can 
reflect underlying processes, such as establishment, growth, 
competition, reproduction, senescence and mortality, which are 
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affected by environmental factors (Sterner et al. 1986). 

The Chittagong Hill Tracts (CHT) region is the major hill 
range of Bangladesh covering about 10% of the total land area of 
the country. It consists of diverse landscapes from hills to water 
bodies and is rich in flora and fauna. Phyto-geographically, it 
shows the admixture of Indo-Chinese floristic elements. The 
forests of this region may be broadly classified into tropical ev¬ 
ergreen, semi evergreen and deciduous types. Commercial tree 
plantations, illegal logging, dam mega-projects, and forced dis¬ 
placement are responsible for the accelerated destruction of the 
precious ecosystems and biodiversity of this region. Rubber, teak 
and eucalyptus monocultures for export have provoked negative 
ecological effects and converted the forest into degraded forest 
land. Despite having only about 10% of the country’s land area, 
the CHT possess one-third of the flowering plant species of the 
country. The people of CHT were abruptly deprived of the tradi¬ 
tional community ownership of lands, which they used for their 
homesteads, jhum cultivation, extraction of forest resources, 
hunting and gathering, etc. 

The government implemented reforestation with Tectona 
grandis L. f. in the degraded areas of CHT. Natural patches have 
grown up in some places of this region where government could 
not take reforestation. Our concern was about the present eco¬ 
logical condition and structure of such a natural patch of forest in 
Kaptai national park of CHT. Therefore, the objectives of our 
study were (1) to elucidate the woody floristic composition and 
species diversity of tropical wet evergreen forest in Kaptai na¬ 
tional park of CHT, and (2) to quantify spatial distribution and 
association of trees. 


Japan for the forest stratification study (Hozumi 1975). There¬ 
fore, our sampling plot was established in the middle portion of 
the mountain as a representative of the forest and had an area of 
900 m 2 (30 m x 30 m) that was divided into 36 quadrats (5 m x 5 
m). The plot, which faced south with a slope of 31.2°, was 77.11 
m above sea level. All woody plants were numbered and identi¬ 
fied to species according to the nomenclature of Das and Alam 
(2001). The measurements obtained for the area were: tree height 
( H) and stem diameter at a height of 10% HI 10 of trees height 
(Arm)- 

Data analysis 


Species dominance 

Dominance of a species was defined by the importance value IV 
(%) which was calculated from the equation proposed by Feroz 
et al. (2008): 
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where H i is the number of individuals of the ith species, Clj is 
the basal area at a height of HI 10 of trees belonging to the ith 
species, f, is the number of quadrats in which the ith species 
appeared, and S is the total number of species. 


Materials and methods 

Study site and sampling plot 

The study was conducted on a tropical wet evergreen forest, 
located at Rampahar (mountain), Kaptai national park 
(22°29.991' N and 92°10.722' E) in the Chittagong hill tracts of 
Bangladesh. This mountain is about 500 m high from the sea 
level (Uddin and Hassan 2012) and was declared as a first re¬ 
serve forest within the Chittagong Hill Tracts (CHTs) in 1875 
(Anonymous 1960). The top part of the mountain is like a mossy 
forest with high humidify. However, the middle part of the 
mountain is really the tropical wet evergreen forest. In this part, 
subsoil is yellowish brown. Loamy soil developed on consoli¬ 
dated or semi-consolidated siltstone, sandstone, and clayey soil 
on shale. Soil on the high hill is excessively drained and less than 
1 m depth. The soil pH ranges from 4.5 to 6.0. The mean 
monthly minimum temperature and the mean monthly maximum 
temperature are respectively 24°C in December and 35°C in May. 
Mean annual temperature was 29.6°C. The mean annual rainfall 
is 2,540 mm. The maximum wind velocity recorded is 96.54 
km/h. 

A sampling plot was less than 400 nr on the Yaeyama Islands, 
Japan for the vegetation survey (Niiro et al. 1974; Suzuki 1979; 
Niiro 1981; Miyawaki 1989), and was 400 m 2 on Iriomote Island, 
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Species-area relationship 

The expected number of species S appeared within the number 

of quadrats ( q ) selected at random from the total number of qua¬ 
drats ( Q ) was calculated from the equation proposed by Shino- 
zaki (1963) (cf. Hurlbert 1971): 




1- 


Q-q^ 
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( 2 ) 


where Clj is the number of quadrats in which the ith species 
occurred. The S q -values were obtained for g-values of 1, 2, 4, 8, 
16, 32 and 36. 

Species diversity 

The following two indices of Shannon’s (MacArthur and MacAr- 
thur 1961) H' and Pielou’s (1969) index j' were used to meas¬ 
ure woody species diversity or equitabilify (evenness). 
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where N is the total number of individuals and the unit of H' is 
bit, or the unit of entropy (e.g. Goldman 2005). 
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^ m A m B 
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Spatial distributions of trees 
* 

The unit-size 771 u - 771 method with successive changes in 
quadrat sizes (Iwao 1972) was used to analyze the spatial distri¬ 
bution of trees. Mean density 171 is defined as: 





(5) 


where X • is the number of individuals in the yth quadrat, and 
is the total number of quadrats when the quadrat size is u. 


In Eq. (7), m* AB quadrat-based mean crowding on species A by 
species B, is defined as: 

o 
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where Q is the total number of quadrat existed in species A and 
species B, and n A j and n Bj are the number of individuals belong¬ 
ing to the jth quadrat in species A and species B, respectively. 

Similarly, m* B A , quadrat-based mean crowding on species B 
by species A, is defined as: 


The C[ u -values were 36, 16, 8, 4, and 2, respectively, for the u- 
values of 1 (5 m x 5 m), 2, 4, 8, and 16. On the other hand, mean 
* 

crowding 771 is defined as the mean number of other individu¬ 
als per individual per quadrat (Lloyd 1967) as follows: 
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If the basic component of the spatial distribution is a single indi¬ 
vidual tree, individual trees are considered to be randomly dis- 

* * 
tributed when 771 u = 771 u , aggregately distributed when 771 u > 

* 

771 u , and uniformly distributed when 771 u < 771 u for any quad¬ 
rat size. 


Interspecific spatial association 

On the basis of the concept of mean crowding proposed by Lloyd 
(1967), Iwao (1977) derived the co-index for analyzing spatial 
associations between species. 
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The interspecies quadrat-based mean crowding indicates the 
mean number of individuals of the other species per individual of 
the subject species per quadrat. On the other hand, m * A , quadrat- 
based mean crowding within species A, is defined as: 


Q 
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Similarly, 7 ?j* b , quadrat-based mean crowding within species 
B, is defined as: 

Q 
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The interspecies quadrat-based mean crowding indicates the 
mean number of the other individuals per individual of the sub- 

Q 

ject species per quadrat. The symbol m A (= IX/ / Q ) and m B 
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A and species B, respectively. The value of to changes from the 
maximum of+1.0 for complete positive interspecies association, 
through 0.0 for independent association, to the minimum of -1.0 
for complete negative association. 

Dendrogram 

The dendrogram for analyzing the degree of spatial association 
of species was constructed following Mountford’s (1962) me¬ 
thod, or unweighted paired-group method using arithmetic aver¬ 
ages (Sneath and sokal 1973). 

Results and discussion 

Floristic composition 

We recorded 25 families, 37 genera, 40 species and 1771 woody 
individuals in the sampling plot. Out of 25 families, 19 families 
(76%) consisted of a single species. Euphorbiaceae and 
Moraceae were the most species-rich families, each of which had 
four species. Out of 37 genera, 34 genera (92%) consisted of a 
single species. Castanopsis, Ficus and Terminalia were the most 
specious genera; each of which contained two species. Of the 40 
species, 10 species (25%) were represented by single individuals. 
Bursera serrata Wall ex Colebr. was the most dominant species 
in terns of the highest importance value (13%) in the total stand 
(Table 1). Therefore, Kaptai national park in Chittagong Hill 
Tracks is dominated by B. setrata. However, this area is highly 
populated by Streblus asper Lour, with the highest number of 
regeneration. Trema orientalis (L.) B1 was typically a light de¬ 
manding species as it appeared only in the top canopy with the 
seventh highest importance value, but had no regeneration (Table 
1). This phenomenon may indicate shade-intolerance for T. ori¬ 
entalis, i.e. early successional species (pioneer species), and also 
indicate a sign of threatening for the present forest. In addition, 
some other species (approximately 27%) with only a few indi¬ 
viduals (one to three) may be only shade species or light de¬ 
manding species which may also be in danger for their surviv¬ 
ability. 

Species-area relationship 

As shown in Fig. 1, the expected number of species increased 
with an increase in the number of quadrats. The relationships of 
the expected number of species S q to the number of quadrats q in 
the total stand were accurately approximated using the following 
equation (Ogawa 1980; cf. Hagihara 1995): 

1 1 1 

— = — r +- (8) 

S, cq d S max 

where c and d are coefficients, and )S max is the expected maxi¬ 
mum number of species. 

The expected maximum number of species was 140 in the to¬ 
tal stand. This result indicates many species may newly invade 
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the forest in the near future as the recorded total number of spe¬ 
cies is 40. Hence, this type of changing species composition, i.e. 
appearance of new species (Fig. 1) and disappearance of species 
(approximately 30% species having only one to two individuals 
in the total stand) (Table 1) in the present forest may express 
instability of floristic composition. 



Fig. 1: Species-area relationship (area of one quadrat: 5 m X 5 m). The 
curve is given by Eq. (8), where c = 9.13, d = 0.50 and S, mx = 144 ( R 2 = 
0.99) for the total stand. 

Species diversity 

The values of Shannon's index H' and Pielou's index J' 
(evenness) were respectively 3.36 bit and 0.63 (Table 2). The 
values of H' and J' in the present forest are higher than the 
values of those reported in a subtropical evergreen broadleaf 
forest in Guangzhou, South China (Wu et al. 2010), whereas 
species richness and expected maximum number of species 5 max 
are lower in the fonner forest than those in the latter forest. It can 
be concluded that the tropical wet evergreen forest in Chittagong 
hill tracts, Bangladesh and the subtropical evergreen broadleaf 
forest in Guangzhou, South China seem to be most likely unequi¬ 
librium in floristic composition with moderately high species 
diversity, because the values of S max are higher than the species 
richnesses for both of the forests. The subtropical evergreen 
broadleaf forest in Okinawa, Japan is stable with very high spe¬ 
cies diversity, because the expected maximum number of species 
(62) is very close to the total numbers of species (60) (Table 2). 

Nevertheless, the values of H' increased rapidly up to 8 (200 
m 2 ) quadrats and then the values tended to be became stable 
despite an increase of area, whereas the values of J' gradually 
decreased with increasing area (Fig. 2). This result suggests that 
the sample area of 200 m 2 in this forest would be sufficient for 
measuring H' . On the other hand, the trend of J' may indicate 
that the rate of equality of individuals among the different spe¬ 
cies decreased with increasing area. 
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Table 1: 

List of 50% woody species of the total stand in order of species rank determined by importance value IV (dominant species) 


Spp. rank 

Scientific name 

Family name 

No. 

of individual 

mi 


Abbreviation 

i 

Bursera serrata Wall ex Colebr. 

Burseraceae 


198 

12.84 


Bs 

2 

Streblus asper Lour. 

Moraceae 


352 

10.55 


Sa 

3 

Castanea indica Roxb. 

Fagaceae 


315 

9.78 


Ci 

4 

Microcos paniculata L. ex W & A 

Tiliaceae 


284 

9.66 


Mp 

5 

Cassia siamea Lam. 

Caesalpiniaceae 


144 

9.18 


Cs 

6 

Maesa ramentacea A. DC 

Myrsinaceae 


190 

7.91 


Mr 

7 

Trema orientalis (L.) B1 

Ulmaceae 


1 

7.45 


To 

8 

Tectona grandis L. f. 

Verbenaceae 


11 

5.21 


Tg 

9 

Meliosma simplicifolia (Roxb.) Walp 

Sabiaceae 


52 

4.55 


Ms 

10 

Lagerstroemia speciosa (L.) Pers. 

Lythraceae 


8 

4.19 


Ls 

11 

Pterospermum acerifolium (L.) Willd 

Sterculiaceae 


74 

3.79 


Pa 

12 

Terminalia crenulata Roth. 

Combretaceae 


4 

3.26 


Tc 

13 

Anogeissus acuminate (Roxb.) Wall, ex Guill. & Perr. 

Combretaceae 


13 

1.33 


Ac 

14 

Vitex glabrata R.Br. 

Verbenaceae 


18 

1.12 


Vg 

15 

Baccaurea ramiflora Lour. 

Euphorbiaceae 


29 

1.09 


Br 

16 

Syzygium fruticosum (Roxb.) DC 

Myrtaceae 


8 

0.95 


Sf 

17 

Toona ciliate M. J. Roem. 

Meliaceae 


11 

0.90 


Tc 

18 

Dioscorea deltoidea Wall, ex Griseb. 

Dioscoreaceae 


12 

0.77 


Dd 

19 

Caesalpinia bonduc (L.) Roxb. 

Caesalpiniaceae 


1 

0.71 


Cb 

20 

Phyllanthus emblica L. 

Euphorbiaceae 


3 

0.38 


Pe 

S, the number of species per sample area 







Table 2: 

ics 

Comparison in species diversity, expected maximum number of species, tree density and basal area among the same type of forests in the trop- 

Sample Area Tree Size Basal Area 

Study Site 7 . , 

(ha) (m) (m 2 h ) 

No. of trees 
(Stem-ha* 1 ) 

S 

Smax 

H' 

J' 

Source 

Okinawa, Japan 0.08 H>0.\ 50.3 

62452 

60 

62 

4.83 

0.82 

Feroz et al. 2006 

Guangzhou, South China 0.12 H> 0.0 35.0 

119836 

75 

127 

3.01 

0.48 

Wu et al. 2010 

Kaptai National Park, BD 0.09 i/>0.1 39.0 

19677 

40 

144 

3.37 

0.63 

Present Study 



Cumilative quadart size 


Fig. 2: Relationships of H' and j' to area. Open circle, If ; closed 
circle J'. The values of H' and ,/' arc the mean of ten samples, which 
were randomly selected. 


Spatial distribution of trees 

The spatial distribution patterns of trees for the total stand and 
top five dominant species are shown in Fig. 3 based on unit-size 
* 

m -m relation (Iwao 1972). The distribution pattern for the 

* 

total stand was completely random, because the M — 7W rela¬ 
tion lines did not differ from the Poisson line. Similarly, Casta- 
nea indica and Cassia siamea were randomly distributed. On the 
other hand, Bursera serrata, Streblus asper and Microcos pani- 
culata showed a different trend in spatial distribution, i.e. at the 
beginning there seemed to exist an aggregate distribution and 
later tended to be random distribution with increasing quadrat 
size. Therefore, dominant species showed random distribution 
for large area. Thus, they may maintain total stand’s random 
distribution. Such type of distribution pattern is probably general 
for the tropical forest. 
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Mean density ra u 

Fig. 3: Relationships between mean crowding m a and mean density m u with successive changes of quadrat sizes for spatial distribution of total stand and 

five dominant species 


Spatial association between species 

The spatial association between species was analyzed by the co¬ 
index (Iwao 1977) for 25% of the total number of species (i.e. 10 
species) in terms of their number of individuals arranged in de¬ 
scending order (Table 3). The strongest positive interspecific 
association occurred between Streblus asper Lour, and Castanea 
indica Roxb. (co =0.51). Microcos paniculata L. ex W & A and 
Meliosma simplicifolia (Roxb.) Walp (co =0.39), Streblus asper 
and Bursera serrata Wall ex Colebr. (co =0.31), and Streblus 
asper and Cassia siamea Lam. (co = 0.28) showed a modest as¬ 
sociation). A total of 58% species associations was completely 
negative i.e. co = -1. As a whole, most species were spatially 
associated weakly with each other. 

The cluster analysis of interspecific spatial pattern of associa¬ 
tion between species produced three major clusters (A, B and C), 
though the degree of spatial association or the degree of overlap¬ 
ping within these clusters was a little less than 0.0 (Fig. 4), i.e. 
clusters A, B and C were slightly exclusive with each other. 
However, almost species in the clusters were modestly associated 
with each other. Cluster A contains four species ( Streblus asper, 
Castanea indica, Bursera serrata and Cassia siamea), cluster B 
contains three species ( Maesa ramentacea, Pterospermum aceri- 
folium and Baccaurea sapida) and cluster C contains three spe- 
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cies ( Microcos paniculata, Meliosma simplicifolia and Vitex 
glabrata). 


Table 3: Spatial association of ten major species in terms of their large 
number of individuals analyzed by the value of co-index 


Species Sa 

Ci 

Mp 

Bs 

Mr 

Cs 

Pa 

Ms 

Bsa 

Sa 0.51 

-0.26 

0.31 

0.09 

0.28 

0.02 

-0.11 

-0.86 

-0.03 

Ci 

-0.18 

0.22 

0.12 

0.07 

-0.11 

-0.04 

-0.48 

-0.58 

Mp 


-0.16 

-0.19 

-0.05 

-0.23 

0.39 

-0.47 

0.15 

Bs 



0.01 

0.15 

-0.16 

-0.26 

0.07 

0.05 

Mr 




-0.17 

0 .18 

-0.30 

0.11 

-0.62 

Cs 





-0.16 

-0.16 

-0.07 

-0.08 

Pa 






-0.08 

0.13 

-0.73 

Ms 







-0.21 

0.24 

Bsa 








0.10 


Note: species abbreviation: Sa — Streblus asper, Ci— Castanea indica, Mp— 
Microcos paniculata, BS — Bursera serrata, Mr— Maesa ramentacea, Cs — 
Cassia siamea. Pa— Pterospermum acerifolium, Ms —Meliosma simplicifolia, 
Bsa — Baccaurea sapida 

The result of cluster analysis showed a general pattern of spe¬ 
cies association, in which species pairs were spatially independ¬ 
ent at all or most small clusters. Species of cluster B and C have 
similar habitat preference, where different species existed, 
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whereas those of cluster A have specific preference, where a 
moderately strong association of species exists. It is likely that 
stands of species from all clusters are mosaics of complete habi¬ 
tat and pioneer habitat. As a whole, it is apparent that all patches 
in this forest community have similar habitat and regeneration 
niches, which could be a phenomenon for a young growth forests. 
Thus, the existence of habitat and regeneration niche may be 
important factor in the maintenance of diversity in this forest. 



Species 

Fig. 4: Cluster diagram showing species spatial association (see Table 1 
for species abbreviation). Sa — Streblus asper, Ci —Castanea indica, Mp— 
Microcos paniculata, BS — Bursera serrata, Mr— Maesa ramentacea, Cs 
— Cassia siamea, Pa— Pterospermum acerifolium, Ms —Meliosma 
simplicifolia, Bsa — Baccaurea sapida 

Tree density and basal area 

Tree density and basal area for D 0 .ih in the present forest were 
19,677 ha' 1 and 39.0 m 2 -ha' 1 , respectively, at the stand level (H> 
0.10 m) (Table 2). These results are lower than the values of 
62,452 ha' 1 and 50.3 nr-ha" 1 (H > 0.10 m) in a subtropical ever¬ 
green broadleaf forest in Okinawa, Japan (Feroz et al. 2006). 
However, the tree density in the tropical wet evergreen forest in 
CHTs, Bangladesh is much lower than that of 119,836-ha" 1 (H > 
0.0 m) in a subtropical evergreen broadleaf forest in Guangzhou, 
South China, whereas the basal area in the fonnaer forest is 
slightly higher than that of 35.0 nr-ha' 1 (H > 0.0 m) in latter 
forest (Wu et al. 2010). The above results suggest that the pre¬ 
sent forest is lowly populated with a little high value for basal 
area, whereas the subtropical evergreen broadleaf forest in 
Guangzhou, South China is highly populated with low value for 
basal area. On the other hand, the subtropical evergreen broad¬ 
leaf forest in Okinawa is well populated with a high value for 
basal area. 
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